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ABSTRACT 

Using simple infrared color selection, the Two Micron All Sky Survey (2MASS) has found a large number of 
red, previously unidentified, radio-quiet quasi-stellar objects (QSOs). Although missed by UV/optical surveys, 
the 2MASS QSOs have Xj-band luminosities that are comparable to "classical" QSOs. This suggests the possible 
discovery of a previously predicted large population of dust-obscured radio-quiet QSOs. We present the results 
of an imaging survey of 29 2MASS QSOs observed with the Wide-Field Planetary Camera 2 onboard the Hubble 
Space Telescope, /-band images, which benefit from the relative faintness of the nuclei at optical wavelengths, 
are used to characterize the host galaxies, measure the nuclear contribution to the total observed /-band emission, 
and to survey the surrounding environments. The 2MASS QSOs are found to lie in galaxies with a variety of 
morphologies, luminosities, and dynamical states, not unlike those hosting radio-quiet PG QSOs. Our analysis 
suggests that the extraordinary red colors of the 2MASS QSOs are caused by extinction of an otherwise typical 
QSO spectrum due to dust near the nucleus. 

Subject headings: quasars: general — galaxies: active — galaxies: interactions — dust, extinction 



1. INTRODUCTION 

Mostly immune from absorption processes and obscuration 
by dust, radio surveys are ideal for identifying complete sam- 
ples of quasi-stellar objects (QSOs). However, only ~10% of 
QSOs are "radio-loud" 1 , necessitating the use of other more 
vulnerable wavelengths. Optical surveys designed to find UV- 
excess objects have long prevailed. Despite their success, 
searches in the optical are fated to miss QSOs that might be 
reddened by dust. 

Webster et al. (1995) studied a large sample of radio-selected 
quasars and found an enormous amount of scatter in their 
Bj-K colors (1 < Bj - K < 8). This was interpreted as the 
result of varying amounts of dust within the quasar host galaxy. 
Optically-selected QSOs showed a much smaller color range, 
presumably because the redder objects were not detected. If the 
dust content around radio-loud and radio-quiet QSOs does not 
differ significantly, their results suggested that approximately 
80% of radio-quiet QSOs are missed by optical surveys. 

The prediction of a large, predominantly undetected popu- 
lation of red QSOs is consistent with theoretical explanations 
for the hard x-ray background (Comastri et al. 1995). However, 
some have claimed that the large spread in Bj - K colors seen 
by Webster et al. (1995) is not due entirely to dust obscuration 
(Boyle & di Matteo 1995; Benn et al. 1998). In that case, any 
missing population of QSOs would be much smaller, if even 
significant. 

The discovery of highly reddened and buried QSOs in the 
Infrared Astronomical Satellite (IRAS) database (e.g., Beich- 
man et al. 1986; Low et al. 1988; Cutri et al. 1994; Hines et al. 
1995, 1999a,b) confirmed that at least some radio-quiet QSOs 
are indeed obscured along our line of sight by dust. While IRAS 

1 Radio-loud QSOs are distinguished by their radio to optical flux ratio, R > 
10. This cut-off value, as well as the wavelengths used, vary throughout the 
literature. 



was only able to detect the mid and far-infrared emission of ex- 
tremely reddened QSOs, it was quickly recognized that a wide 
area, near-infrared survey could potentially uncover large num- 
bers of previously unknown QSOs due to less extinction at 2/im 
than in the optical. 

With this key feature in mind, Cutri et al. (2001) are compil- 
ing a catalog of active galactic nuclei (AGN) derived from the 
Two Micron All-Sky Survey (2MASS: Skrutskie et al. 1997) 
Point Source Catalog. Even with a relatively small fraction of 
the sky analysed so far, spectroscopic followup has confirmed 
over 400 previously unidentified emission-line AGN with QSO- 
like luminosities (M Ks < -23 and z < 0.4: Cutri et al. 2001, 
2002). 

The nature of these red 2MASS QSOs is an open ques- 
tion. The ZT-band and [O III]A5007 luminosities of many of the 
newly found objects fall in the same ranges as for UV/optically- 
selected QSOs (Smith et al. 2003). However, nearly all have 
much redder colors (B-Ks) compared with traditional QSOs 
(e.g., PG QSOs: Schmidt & Green 1983), implying that many 
of the 2MASS objects are obscured from our direct view by 
dust. Smith et al. (2002) find that at least 10% of the 2MASS 
QSOs are "highly-polarized" (P > 3%); a significantly higher 
fraction than both the PG (Berriman et al. 1990) and broad ab- 
sorption line QSOs (Schmidt & Hines 1999). Indeed, only the 
/RAS-selected hyperluminous infrared galaxies (HIGs) show 
higher polarization. In addition, several of the highly polar- 
ized 2MASS QSOs exhibit broad polarized emission lines even 
though they are dominated by narrow emission lines in total 
flux (Schmidt et al. 2002; Smith et al. 2003). This is very remi- 
niscent of the simple orientation-dependent unified schemes for 
Seyfert galaxies (Antonucci 1993; Tran 1995) and HIGs (Hines 
& Wills 1993; Hines et al. 1995, 1999a,b; Goodrich et al. 1996; 
Tran 2001). In these schemes, a dusty torus obscures our direct 
view but light from the nucleus escapes through the open poles 
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to be scattered (thus polarized) into our line of sight; type 2 ob- 
jects appear different from type 1 simply because their tori are 
more highly inclined. 

However, there is still an ambiguity between simple orien- 
tation schemes and those that involve significant dust cover 
from debris left by strong galaxy-galaxy interactions. It has 
been suggested that ultraluminous infrared galaxies (ULIRGs), 
//?A5-selected systems with the bulk of their energies emerg- 
ing in the far-infrared and with bolometric luminosities of 
10 12 -10 13 L 0, can evolve into classical QSOs once the gas 
and dust around their nuclei are either consumed as fuel by a 
powerful starburst or a nascent AGN, or blown away in pow- 
erful winds (Sanders et al. 1988; Heckman, Armus & Miley 
1990). In this scenario, some or all of the 2MASS QSOs might 
be aging ULIRGs, complete with post-starburst stellar popula- 
tions, fading tidal remnants, disturbed morphologies, and per- 
haps closely spaced multiple nuclei. 

The 2MASS-selected QSOs might therefore be (1) relatively 
normal QSOs seen from near the plane of an obscuring torus, 
(2) fairly young, evolving post-ULIRG systems wherein most 
of the molecular gas and dust has been removed from the cir- 
cumnuclear environment, but that still retain enough material 
to make them "invisible" to classical UV or visual searches for 
AGN, or (3) an entirely new aspect of the overall phenomenon 
of nuclear activity. As part of a multi-wavelength investiga- 
tion of the 2MASS AGN, we have performed a snapshot imag- 
ing survey with the Wide-Field Planetary Camera 2 (WFPC2) 
aboard the Hubble Space Telescope (HST) to investigate these 
possibilities. The images have enabled us to characterize the 
host galaxies, estimate the relative contributions to the emitted 
optical flux from the host galaxy and the active nucleus, and 
survey the immediate environments of these objects. 

2. DATA 

2.1. The 2MASS QSO Sample 

Candidate AGN were selected from 2MASS according to 
color (J-Ks > 2), galactic latitude (\b\ > 30°), and detection 
(complete to K s < 15.0) in each of the three 2MASS bands, 
JHK S (Cutri et al. 2001, 2002). These criteria efficiently fil- 
ter almost all known AGN and all but the coolest stars (Beich- 
man et al. 1998). Approximately 80% of these candidates were 
confirmed as AGN via follow-up optical spectroscopy that also 
determined their redshifts and spectral type. Of these AGN, 
most have absolute K s magnitudes that easily fall within the lu- 
minosity range of UV/optically-selected QSOs (Figure 1). A 
representative subset (90) of this red, radio-quiet QSO (M Ks < 
-23) catalogue was surveyed for optical broad-band polariza- 
tion (Smith et al. 2002). Twenty-nine of these objects comprise 
the sample addressed by this paper. They range in redshift from 
z = 0.136 to z = 0.596 with a median value of 0.213. The sample 
members and their previously determined properties are pro- 
vided in Table 1 . For brevity, the designation of objects from 
the 2MASS Incremental Release Point Source Catalog has been 
shortened to 2Mhhmmss±ddmm in the text. 

2.2. Observations 

A WFPC2 snapshot survey was conducted between the UT 
dates of 2 July 2001 and 14 October 2001 during Cycle 10 of 
the HST observing schedule. Twenty-eight of fifty-four poten- 
tial targets were observed during that period. A final twenty- 
ninth object (2M0 10607+2603) was added on 10 June 2002. 



The QSOs were centered on the planetary camera (PC) and ob- 
served with a single filter, F814W. Each observation consisted 
of two 400 second integrations in CR-SPLIT mode to allow for 
the removal of cosmic rays. In order to maximize integration 
time, no dithering was attempted. The survey was designed to 
balance the scientific requirements of the project with the sim- 
plicity needed for successful snapshot scheduling. WFPC2 was 
selected to facilitate comparisons with previous studies, and the 
PC, with 0.046" pixels and a 36.8" field of view, was cho- 
sen over the wide field cameras in order to better sample the 
point spread function (PSF). As a result, the light contribution 
of AGN can be subtracted more precisely, but at the expense 
of sensitivity to extremely low surface brightness features. The 
F814W filter, which corresponds closely to the Cousins /-band 
(hereafter referred to simply as I), has the dual advantage of be- 
ing sensitive to red light from older stellar populations, while 
still benefitting from the relative faintness of the nuclei at opti- 
cal wavelengths. 




-24 -26 -28 -30 

FIG. 1. — B-K s color vs. absolute K s magnitude. The 2MASS QSOs are 
much redder than the UV-excess PG QSOs but have comparable near-infrared 
luminosities (Neugebauer et al. 1987). The circles are a representative in- 
ventory of the red QSOs discovered by 2MASS. Those investigated in this 
paper are indicated by ®. The PG QSOs plotted here are the 75 (out of 1 14) 
with Mk < —23 and z < 0.6. The reddening vector shown in the upper right- 
hand corner represents one magnitude of extinction in the V-band (Ag = 1.324, 
A K = 0.1 12: Rieke & Lebofsky 1985). 

2.3. Data Reduction 

The WFPC2 images were calibrated using the STSDAS rou- 
tine calwp2 that is incorporated into the STScI data reduction 
pipeline and makes use of the best reference files available. 
Subsequently, a flat-field correction was applied as described in 
the WFPC2 Instrument Science Report 2001-07 (Karkoshka & 
Biretta 2001). The two images for each object were combined 
and cosmic rays were removed using the IDL routine cr_reject 
that is nearly identical to crrej in IRAF and is available in the 
IDL Astronomy User's Library. No compensation was made 
for charge transfer efficiency (CTE) problems known to affect 
WFPC2, resulting in an underestimation of the illumination in 
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Table 1 

Sample of 29 2MASS-selected red QSOs 
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"Redshifts from Cutri et al. (2001). 

*Near-infrared apparent magnitudes and colors from the 2MASS Point Source Catalog. 
c Calculated (Eq. 1) from Column 5 values (Hq = 75 km s _1 Mpc~' and qo = 0.5). 
''Spectral classifications from Cutri et al. (2001) and Smith et al. (2003). 

^Following recent observations, this object no longer meets the ]-K s > 2 selection criterion. It has been included in the analysis of the sample. 



the center of the chip at the few percent level (Whitmore & 
Heyer 1997; Biretta et al. 2002). The dominant source of error 
in the magnitudes we derive is not the CTE compensation, but 
uncertainty in the PSF-subtraction. Sky subtraction was carried 
out using custom IDL software. 

2.4. PSF -Subtraction 

Factors affecting the PSF include the filter response, scat- 
tering due to surface irregularities in the optics such as dust, 
the shape of an object's spectral energy distribution (SED), and 
especially slight changes in the position of the secondary mir- 
ror, or "breathing". While all of these may vary slightly be- 
tween observations, their combined effect is minimized by cre- 
ating a PSF directly from observed point sources. Using the 
sample QSOs themselves for this purpose is preferable to field 
stars for two reasons. The PSF is somewhat dependant on a 
source's position on the chip, and a star's color significantly 
affects whether or not it is suitable for modeling a QSO PSF 
(Bahcall et al. 1997). As this observing program was designed 
to be a single filter snaphot survey, colors were not available and 
valuable orbits were not dedicated to observing stars centrally 



located on the chip. 

In order to identify the most "point source like" QSOs in our 
sample, an artificial PSF created by the software package Tiny 
TIM (Krist & Hook 2001) was scaled and subtracted from each 
image where a host galaxy was not readily apparent. Direct 
comparison of the resulting profiles revealed extended emis- 
sion around six of the nine candidates. Two of the remaining 
three were not optimal for the purpose of constructing a PSF. 
2M 170003+21 18 was relatively faint requiring that it and its in- 
herent noise be scaled significantly to combine it with brighter 
objects. A saturated star in the field of 2M023430+2438 was 
a source of contamination that could not be removed reliably. 
Therefore, our PSF was created solely from 2M222202+1952. 

The light contribution of the AGN can be bounded by consid- 
ering two extreme cases. The upper limit to the AGN emission 
corresponds to the PSF amplitude which results in the subtrac- 
tion residuals summing to zero within a specified aperture (de- 
fined here as the full width at half maximum (FWHM) of the 
PSF). With the constraint that the PSF amplitude must be suffi- 
ciently large to eliminate the presence of diffraction spikes and 
an Airy pattern, the lower limit corresponds to the removal of 
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the central peak without introducing any negative residuals. 

Between these bounds, the PSF amplitude which minimizes 
the variance of the residuals can be considered the best deter- 
mination of the light contribution of the AGN. This approach 
yields generous errors. Finally, for the purpose of estimating a 
3cr flux uncertainty due to PSF-subtraction, the smaller of the 
differences between the accepted value and the upper and lower 
limits was adopted. 



2.5. QDOTULIRGs 

Before presenting the results from the analysis of the 2MASS 
QSOs, we briefly describe a second sample of objects inves- 
tigated in this paper. For reasons outlined in §1 we wished 
to compare properties of the 2MASS QSOs to a sample of 
ULIRGs. Identification of companion objects surrounding the 
latter are not currently available in the literature, requiring us to 
perform our own analysis of HST archival data. A substantial 
number of ULIRGs have been observed with WFPC2/F814W 
(Borne, Bushouse, Lucas, & Colina 2000) and are publicly 
available in the archive. These images were taken in exactly 
the same manner as the 2MASS QSO sample, with the excep- 
tion of being centered on WF3 instead of the PC. Cross-listing 
those objects with the QDOT 2 ULIRG sample and imposing the 
same sky position restriction used to select the 2MASS QSOs 
(\b\ > 30°) results in a homogeneous sample of 30 high-latitude 
ULIRGs. These objects and their redshifts are presented in 
§3.5. Precisely the same steps described in §2.3 were applied 
to these observations. 



3. RESULTS 

Images of all 29 2MASS QSOs prior to and following PSF- 
subtraction are shown in Figure 2. The entire sample is dis- 
played at the same angular and physical scales in the first and 
second columns, respectively. The second and third columns 
show the objects before and after subtraction of the light contri- 
bution from the AGN. In some cases this difference is dramatic; 
in others there is very little visual change. The position angle of 
polarization (taken from Smith et al. (2002, 2003)) for those ob- 
jects with measured P jap > 4 is denoted by the double-headed 
arrow in the upper righthand corner in column two. Isophotes 
of the residual images (and the masks used in calculating them) 
are shown in the last column. In some cases, the isophotes more 
clearly reveal underlying structure, highlighting features such 
as multiple nuclei and isophotal twisting caused by interaction 
or the presence of spiral arms. The solid contours indicate flux 
levels that are distinguishable from image noise at the 3cr con- 
fidence level, as opposed to the dotted contours. The dashed 
line identifies the isophote used to determine the eccentricity 
and position angle of the host galaxy. Measured properties of 
the host galaxies are listed in Table 2 and are discussed be- 
low. Additional notes on individual objects can be found in 
Appendix A. 

2 QDOT is a redshift survey of 2,387 IRAS galaxies brighter than the IRAS PSC 
60/im completeness limit (Lawrence et al. 1999, and references therein). In the 
context of this paper, it is only used as means for constructing a ULIRG sample 
with uniform selection criteria. 



3.1. Luminosities 

The PSF-subtraction revealed the relative flux contributions 
of the AGN and host galaxies within 10" radius apertures 3 . Ab- 
solute magnitudes were then calculated using the standard for- 
mula (Peebles 1993) 

m-M = 42.38 + 51og[z/7?] + 1.086(1 -q )z+K (1) 

and the prescription found in the WFPC2 Instrument Handbook 
(Biretta et al. 2002) 

m, = -2.5 log [2.508 x 10 -18 (DN/Exp)] -22.31 (2) 

where h = // /100, K is the redshift-dependant K-correction, 
DN is the number of counts in the F814W filter, and Exp is the 
exposure time. Throughout this paper, values for the cosmo- 
logical constants qo and Ho are taken to be 0.5 and 75 km s" 1 
Mpc" 1 , respectively. Although transformation between the HST 
system and standard magnitudes generally depends strongly on 
the shape of an object's spectrum, the effect is minimal for the 
F814W filter (^0.05 magnitudes). Equation 2 yields results 
nearly identical to those obtained by Holtzman et al. (1995). 

K-corrections for the host galaxies, ranging between -0.21 
and -1.01 magnitudes, were measured from spiral and ellip- 
tical SED templates. For those hosts with ambiguous mor- 
phologies, an average of the two templates was used. For 
the AGN, K-corrections were taken from the Sloan Digital 
Sky Survey (SDSS) composite QSO spectrum (Vanden Berk 
et al. 2001). Of course, the 2MASS QSOs are much red- 
der than UV/optically selected QSOs and their physical na- 
ture has yet to be adequately determined. Extrapolating their 
SED from a "typical" QSO spectrum in order to calculate K- 
corrections is not necessarily correct. Fortunately, for the red- 
shift range of this sample, these corrections are relatively small 
(0. 14 <K< 0.28, (K) = 0.1 8). 

Absolute magnitudes for the host galaxies and AGN are 
given in Table 2. For the latter, both absolute rest frame /-band 
and apparent (without K-corrections) F814W values are listed. 
In this cosmology, an L* galaxy has an /-band absolute magni- 
tude of Mi ~ -22.3 (Norberg et al. 2002; Blanton et al. 2003). 
The 2MASS hosts range between L~ 0.3L* and L- 4L* with 
three quarters of them being brighter than L*. The range in lu- 
minosity of the nuclei extends over 7 magnitudes, with only ten 
of the 29 objects satisfying the traditional M B < -21.5 + 5 \ogh 
criterion 4 for QSO classification (Schmidt & Green 1983). Sec- 
tion 5 presents evidence that AGN luminosity is correlated with 
the Inuc-Ks color index, suggestive of a reddening effect. 

3.2. Surface Brightness Profiles 

Isophotes were fit to the host galaxy images after the point 
source was subtracted and nearby stars/galaxies were masked. 
The methodology used is described by Buote & Canizares 
(1994) and references therein. Surface brightness profiles were 
then constructed using these isophotes. The average magnitude 
per square arc second was calculated as a function of radius 
along the semi-major axis for the regions between successive 
elliptical annuli. Figure 3 shows surface brightness profiles for 
each object, plotted against r as well as r 1 / 4 . Linear profiles on 
the left reflect an exponentially decreasing surface brightness (a 
disk). On the right, linear profiles indicate a de Vaucouleurs r 1 / 4 
law that is characteristic of spiral bulges and elliptical galaxies. 

3 Foreign sources such as stars, background galaxies and companions were 
masked. 

4 Extrapolating to /-band yields M/ < -22.3. 
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FIG. 2. — WFPC2 F814W images of 2MASS QSOs before (first and second columns) and after (third and fourth columns) PSF-subtraction. The position angle 
of polarization for those objects with measured P/ap > 4 is denoted by the double-headed arrow in the upper righthand corner in column two. Isophotes of the 
residual images (and the masks used in calculating them) are shown in the last column. The solid contours indicate flux levels that are distinguishable from image 
noise at the 3(7 confidence level, as opposed to the dotted contours. The dashed line identifies the isophote used to determine the eccentricity and position angle of 
the host galaxy. 
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FIG. 3. — 2MASS QSO host galaxy surface brightness profiles corresponding to the PSF-subtracted images shown in column three of Figure 2. A detailed account 
of these profiles is given in §3.2. 
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The profiles of the three objects with no obviously discernible 
host galaxies were noise-limited and are not included. Addi- 
tionally, it should be noted that the profile shown for the last 
object, 2M234449+1221, may be misleading. The inner points 
are dominated by the immediate region surrounding the QSO, 
while at larger distances the surface brightness reflects the other 
feature clearly seen in Figure 2. Whether this is a double nu- 
cleus within the same host galaxy or a separate companion is 
unclear. By chance this results in a fairly smooth surface bright- 
ness profile. 

3.3. Morphologies 

In three cases, a host galaxy could not be visually distin- 
guished from either the image noise or the PSF-subtraction 
residuals. Considering the survey's short exposure times, the 
PC's relative insensitivity to low surface brightness features, 
and the comparatively high redshifts (z = 0.310, 0.366 and 
0.596) of these objects, it is reasonable to assume that host 
galaxies are present but very faint 5 . In fact, luminosities from 
the residual flux are consistent with the other host galaxies' 
absolute magnitudes. For the remainder of the sample, mor- 
phologies were assigned by visual inspection of both the PSF- 
subtracted images and the surface brightness profiles. In many 
cases, classifying the type of host galaxy is an extremely un- 
certain exercise. Table 2 lists only fairly unambiguous mor- 
phologies, leaving nearly half the sample unclassified. Of the 
26 QSOs with visible hosts, nine reside in disk (spiral) galaxies 
and four are found in what are well described as ellipticals. 

The hosts classifed as spirals preferentially lie at low red- 
shifts (Figure 4). A Kolmogorov-Smirnov (KS) test yields only 
a 0.31% probability that the redshifts of the spirals and remain- 
ing galaxies are drawn from the same parent population. This 
suggests that some, or all, of the unclassified host galaxies may 
have spiral disks that are less obvious due to their smaller an- 
gular extent and the shallow exposures of the shapshot survey. 
This is a testable hypothesis, at least in a coarse sense. Fig- 
ure 5a shows the normalized distributions of axial ratios for 
12,045 spiral and 1,022 elliptical galaxies in the Third Refer- 
ence Catalog of Bright Galaxies (RC3: de Vaucouleurs et al. 
1995). The axial ratios of the unclassified 2MASS QSO hosts 
are shown in Figure 5b. There is a ~ 3% and ~ 63% probability 
that the unclassified host axial ratios are drawn from the distri- 
butions of the RC3 elliptical and spiral galaxies, respectively, 
as indicated by KS tests. Both of these results are inconclusive, 
except to say that it is more likely that the unclassified hosts 
are all spirals than that they are all ellipticals. Their nature can 
be investigated more generally with a Monte Carlo approach. 
A total of N spiral axial ratios and 12 —N elliptical axial ratios 
were selected at random using the RC3 distributions as parent 
populations. These 16 values were then binned to the same 
degree as Figure 5. Repeating this exercise 10,000 times re- 
sulted in an average distribution for N spiral galaxies as well 
as a Gaussian error for each bin. Computing \ 2 indicated how 
well the unclassified 2MASS QSO hosts can be described as a 
sample of N spiral galaxies and remaining ellipticals. This test 
was not very sensitive to the precise choice of N, but it showed 
that the unclassified hosts can be characterized equally well as 
being purely spiral galaxies or evenly divided between the two 
morphologies. 

The preceding analysis begs the question of whether or not 

5 The 800 second observations reached a limiting magnitude of approximately 
Mi < 24. 



the 2MASS QSOs are exclusively found in spiral galaxies. For 
this to be the case, the four hosts classified as ellipticals would 
in fact have to be bulge-dominated spiral galaxies for which the 
disk was simply missed. However, two of these lie at fairly low 
redshifts (z = 0.140 and z = 0.199) where disks would likely 
be detected if present. Also, two-thirds of the 2MASS QSOs 
with unclassified hosts also exhibit evidence of interaction in- 
cluding tidal debris, obviously merging galaxies and general 
asymmetry. This is true of only two other hosts, one spiral 
and one elliptical. Thus, what is seen clearly at low redshift 
as two closely separated galaxies, might instead appear to be 
a nearly edge-on spiral galaxy when seen at a higher redshift 
(e.g., 2M034857+2155 and 2M010607+2603). The results of 
this snapshot survey seem to indicate that 2MASS QSOs are 
found in both spiral and elliptical host galaxies (possibly more 
often in the former) exhibiting a range of dynamical activity. 
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FIG. 4. — Redshift distributions of 2MASS QSOs investigated in this paper. 
The dotted line indicates the subset of values for the host galaxies identified as 
spirals. 

3.4. Inclinations 

The distribution of inclinations for the disk galaxies in the 
sample can potentially reveal selection effects caused by ob- 
scuring dust local to the nucleus or in the host galaxy. Incli- 
nations were calculated while accounting for the intrinsic axial 
ratio, a, of an edge-on galaxy using 



cos; = Vt^' (3) 

where q is the axial ratio b/a and a = 0.2 (Hubble 1926; Holm- 
berg 1958). The inclinations for each host are listed in Table 2 
and their distribution is shown in Figure 6. The mean inclina- 
tion is 67°, with only one inclination greater than 75° and none 
less than 50°. 

Despite our limited sample size (nine), there is an apparent 
disparity between the inclinations of the 2MASS QSO hosts 
and those of UV-selected AGN. With the exception of hard X- 
ray selected samples, spirals hosting AGN preferentially have 
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* Apparent magnitude in HST F814W filter (not K-corrected). 

TJouble-nuclei, merging galaxies or disturbed morphologies such as tidal debris and obvious asymmetry. 
d The sample PSF was created from 2MASSI J222202.2+195231. 
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FIG. 5. — left - (a) Distributions of axial ratios for 12,045 spiral and 1,022 elliptical galaxies in the RC3 Catalog (de Vaucouleurs et al. 1995). (b) Measured axial 
ratios of 2MASS QSO host galaxies with unclassified morphologies. 

FIG. 6. — right - Spiral host galaxly inclinations of 2MASS QSOs and UV-selected AGN (Simcoe et al. 1997), as well as the expectation for randomly oriented 
disks. The area under the curves have been normalized to one. Therefore, the vertical axis represents the fractional number of inclinations in each bin. 
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low inclinations (Simcoe et al. 1997). Maiolino & Rieke (1995) 
found a more even distribution of inclinations in a sample of 
very nearby Seyferts ((d) = 34 Mpc) in which nuclei with low 
apparent brightness were observable. They suggested that a 100 
pc-scale torus coplanar with the galactic disk could be responsi- 
ble for obscuring nuclear activity in highly-inclined disk galax- 
ies. This is consistent with the results of our study. Reddened 
AGN in highly-inclined spirals that are missed by UV/optical 
surveys could be detected in the near-infrared. Furthermore, 
low inclinations would not be found in our sample because 
a less obstructed view of the nucleus would not satisfy the 
2MASS red AGN color criterion. 

A lack of the highest inclinations is not surprising. At i = 90° 
the nucleus would be viewed through the entire plane of the 
galaxy and not likely seen at all. Also, the 2MASS selection 
criteria require objects to be seen in all three wavebands. Ex- 
tremely red "J-band drop-outs" that were excluded from the 
original 2MASS AGN sample, may include objects with host 
galaxies seen nearly edge-on. 

3.5. Apparent Companions 

Previous studies have shown that low-redshift QSOs gener- 
ally lie in moderate galaxy groups as opposed to rich clusters 
or field densities (Hartwick & Schade 1990; Bahcall & Chokshi 
1991; Fisher et al. 1996, and references therein). The environ- 
ment in which this sample of 2MASS-selected QSOs is found 
was characterized by studying the frequency and distribution of 
apparent companion galaxies. This refers to all extended ob- 
jects with small projected separations. Redshift and color infor- 
mation are unavailable and only magnitude considerations can 
be used to eliminate some background galaxies. Clearly, the 
number of apparent companions is greater than or equal to the 
number of true companions that, by definition, lie in close prox- 
imity to the QSOs. The presence of physically unrelated objects 
will also dilute any non-random radial distribution (in the plane 
of the sky) of companion galaxies. However, for the purpose 
of comparative studies, background galaxies do not pose a sys- 
tematic problem as long as the same identification criteria are 
employed for all studies. 



Table 3 

2MASS QSO apparent COMPANIONS 
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projection 
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Table 4 



QDOT ULIRG apparent COMPANIONS 
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-21.0 


8.5 



Apparent companions were selected in the following way. 
The images were marginally smoothed and high-density re- 
gions with surface brightness exceeding the image noise by 
> 3a were flagged. The FWHM was then measured for each 
candidate and compared to that of the previously constructed 
PSF as well as with obvious field stars. Objects with a FWHM 
less than the adopted two pixel cutoff, as well as those pos- 
sessing diffraction spikes in the case of stars with saturated 
cores, were rejected. Finally, an absolute magnitude criterion 
(Mi < -18.3) was imposed assuming the same redshift as the 
QSO. This serves to discriminate against more luminous ob- 
jects at larger, unrelated distances. 
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FIG. 7. — Polarization corrections and 3a error bars for all 29 AGN due to dilution from host galaxy starlight. Dashed lines indicate correction factors of two and 
ten. Each symbol corresponds to an AGN spectral type bin (the E3 refers to the 2/S object). The corrections correlate with AGN type (type 1 through 2). 



The identical procedure was performed on the sample of 30 
ULIRGs drawn from the HST archive. To account for the dif- 
ferent pixel scale of the WF3 camera, these images were first re- 
sampled to match the 2MASS QSO images observed with the 
PC. This allowed us to select the high density surface bright- 
ness regions with the same criteria, preventing any bias that 
might have resulted from different angular pixel sizes. The 
increased sensitivity of the wide field cameras to low surface 
brightness should not present a problem since the magnitude 
cutoff was substantially brighter than the limiting magnitudes 
of the 2MASS QSO images. 

Apparent companions found within 25 h~ l projected kpc of 
the 2MASS QSOs are presented in Table 3. Table 4 provides 
complimentary information for the QDOT ULIRGs in addition 
to their redshifts. Further discussion of these companions is 
deferred to §5 where a third sample will be considered. 

3.6. Polarization 

The optical broadband linear polarization of 70 2MASS 
QSOs (including all 29 objects addressed by this paper) has 
been measured, and a detailed discussion of their properties can 
be found in Smith et al. (2002). More than 10% of 2MASS 
QSOs show high polarization (P > 3%), ranging up to 11%. 
Lower polarization values were generally observed for lower 
luminosity AGN, a trend that can be understood in part by dilu- 
tion of polarized light by unpolarized light from the host galaxy. 

The PSF fitting in §2.4 provides a measure of the polarization 
dilution via starlight. If the linear polarization, P, is measured 
with an uncertainty o>, and the total flux, Ftot, has nuclear 
contribution Fagn with uncertainty <Jf agn , then the nuclear po- 
larization, Pagn, an d uncertainty, <Jp AGN , are given by 



VPagn ~ 



Ftot 
Fagn 



a P + (J F AGN 



Fagn 



1/2 



(4) 
(5) 



Columns 2 and 4 of Table 5 list the measured and derived val- 
ues. No adjustment for statistical bias has been made. 

Figure 7 shows the range of polarization corrections (= 
Ftot /Fagn) with 3cr error bars for the sample of 2MASS QSOs. 
The error- weighted mean correction is 1.46, although 52% 
(15/29) of the corrections are at least a factor of two, and 28% 



(8/29) are ten or greater even with consideration of the mar- 
gin of error. Clearly, the flux contribution from host galaxies 
represents an important correction for the polarization of many 
2MASS QSOs. Without this correction, 17% (5/29) of the sam- 
ple had measured polarizations in excess of 3%, all significant 
at greater than the 4er level. With the correction, 62% (18/29) 

Table 5 

Polarization corrections for starlight dilution 



Object 


pa 


rAGNb 
J JOT 


Pagn 


(2MASSI J) 


(%) 


(%) 


000810.8+135452 


0.79±0.47 


0.01±0.001 


78.83±47.70 


005055.7+293328 


2.47±0.49 


0.18±0.009 


13.42±2.74 


010607.7+260334 


6.03±0.54 


1.00±0.020 


6.03±0.55 


015721.0+171248 


1.44±0.47 


0.40±0.015 


3.64±1.20 


022150.6+132741 


0.39±0.24 


0.08±0.002 


4.60±2.83 


023430.6+243835 


2.57±0.46 


0.90±0.003 


2.85±0.51 


032458.2+174849 


1.23±0.38 


0.76±0.009 


1.62±0.50 


034857.6+125547 


2.23±0.73 


0.12±0.014 


19.36±6.74 


092049.0+190320 


1.15±0.52 


0.23±0.014 


5.07±2.31 


125807.4+232921 


1.22±0.13 


0.87±0.006 


1.40±0.15 


130700.6+233805 


2.45±0.63 


0.07±0.006 


36.10±9.80 


145331.5+135358 


0.81±0.29 


0.07±0.005 


10.80±3.94 


151621.1+225944 


1.02±0.28 


0.55±0.022 


1.85±0.51 


152151.0+225120 


1.02±0.48 


0.07±0.003 


14.18±6.71 


154307.7+193751 


1.33±0.26 


0.49±0.020 


2.72±0.54 


163700.2+222114 


2.34±0.40 


0.35±0.002 


6.70±1.15 


163736.5+254302 


0.42±0.28 


0.25±0.019 


1.65±1.11 


165939.7+183436 


6.30±0.73 


0.56±0.016 


11.21±1.34 


170003.0+211823 


11.11±0.80 


1.00±0.015 


11.11±0.82 


171442.7+260248 


0.86±0.33 


0.57±0.036 


1.50±0.58 


171559.7+280717 


6.08±1.28 


0.56±0.011 


10.81±2.29 


222202.2+195231 


7.19±1.14 


0.99±0.003 


7.26±1.15 


222221.1+195947 


0.95±0.23 


0.76±0.009 


1.25±0.30 


222554.2+195837 


1.38±0.56 


0.01±0.001 


98.82±40.65 


225902.5+124646 


1.03±0.78 


0.49±0.016 


2.09±1.59 


230304.3+162440 


0.24±0.51 


0.02±0.003 


10. 85 ±23. 09 


230442.4+270616 


0.24±0.55 


0.26±0.016 


0.94±2.15 


232745.6+162434 


1.08±0.57 


0.06±0.002 


18.89± 10.00 


234449.5+122143 


1.01 ±0.24 


0.69±0.002 


1.46±0.35 



"Linear polarization measurements (not corrected for statistical bias) from Smith 
et al. (2002). 

^Fractional contribution of AGN to total light (within 10" radius aperture) in 
HST filter F814W. 
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have P > 3% though only 24% (7/29) have both P > 3% and 
P/a P >4. The AGN is likely more obscured in the /?-band (ap- 
proximately where the polarization measurements were made) 
than in the /-band (where the starlight dilution was measured), 
so the correction factors may be even larger than stated. 

The 2MASS QSOs reflect the full range of AGN types from 
type 1 to type 2 (Table 1). Using the spectral classification 
adopted by Cutri et al. (2001), type "1.x" refers to objects 
that exhibit broad Ha emission lines but have no detected 
H(3. Objects that show broad Ha but whose spectral cover- 
age does not include H(i are labeled "Lb". Classification of 
2M232745+1624 is ambiguous. A relationship is evident be- 
tween AGN type and the degree of dilution due to host galaxy 
starlight (Figure 7). Type 1 and 1.2 AGN have the smallest 
Ftot/Facn ratios while type 1.8-2 are dominated by the host 
galaxy. This result is consistent with the assertion by Smith et 
al. (2002) and the findings of Schmidt et al. (2002) and Smith 
et al. (2003) that the low polarization of type 1.8-2 objects is 
largely caused by the contribution of unpolarized starlight. 

4. COMPARISON STUDIES 

The aim of this snapshot survey was to characterize the host 
galaxies and environments in which the red 2MASS QSOs are 
found. Directly comparing these properties to those of previ- 
ously investigated samples offers insight into the nature of the 
2MASS QSOs and allows us to address the questions raised in 
§ 1 . Three samples of ULIRGs were selected from the litera- 
ture, as well as five studies of QSOs and their host galaxies. 
They were primarily chosen on grounds of suitability for useful 
comparison and are, therefore, a representative selection rather 
than an exhaustive one. 

Surace et al. (1998) used HST to study nine ULIRGs, em- 
ploying the same instrumental setup as this 2MASS QSO snap- 
shot survey (WFPC2/PC/F8 14W). Color information allowed 
these authors to investigate a number of issues which are be- 
yond the scope of discussion in this paper. However, they do 
report the incidence of interaction among their sample. 

Farrah et al. (2001) imaged 23 QDOT ULIRGs with HST's 
F606W filter (-V-band) and the PC chip of WFPC2. They 
present observed properties including morphology and interac- 
tion. Only twelve of their objects have \b\ > 30°, a necessary 
condition for fairly comparing their surrounding environoments 
to those of the 2MASS QSOs. For this reason, apparent com- 
panions were identified for a different (but overlapping) sample 
of ULIRGs as was discussed in §2.5. 

Kim, Veilleux, & Sanders (2002) recently obtained ground- 
based optical (/?-band) and near-infrared (/f'-band) imaging of 
the IRAS 1-Jy sample of 1 18 ULIRGs. The analysis of these im- 
ages is presented in Veilleux, Kim, & Sanders (2002). Among 
other things, these authors address morphology and the inci- 
dence of interaction. 

The five QSO samples were constructed with varying selec- 
tion criteria and include significant overlap. To ensure a ho- 
mogeneous data set, only the radio-quiet PG QSOs were con- 
sidered. We briefly describe each study before discussing the 
aggregate PG QSO sample. 

McLeod & McLeod (2001) observed sixteen radio-quiet 
QSOs in the //-band with the Near-Infrared Camera and Multi- 
Object Spectrometer (NICMOS) on HST. Ten of these objects 
were originally part of a sample of 26 high-luminosity QSOs 
(M B < -23. 1 and z < 0.3) selected from the Bright Quasar Sur- 
vey (Schmidt & Green 1983) and studied by McLeod & Rieke 



(1994b). The additional six objects were drawn from a slightly 
larger redshift range (z < 0.4). 

Percival et al. (2001) obtained deep Zf-band ground-based 
imaging for a sample of 14 high luminosity, radio-quiet QSOs. 
The observations were made with the UK Infrared Telescope 
(UKIRT) equipped with the IRCAM3 camera. Their sample 
has total absolute magnitudes of My < -25 and span the redshift 
range 0.26 < z < 0.46. 

Surace, Sanders, & Evans (2001) carried out a ground- 
based optical and near-infrared study of 17 "infrared-excess" 
PG QSOs (all PG QSOs with z < 0.16 and L ir /L BBB 6 > 0.46) in 
order to investigate whether ULIRGs represent the progenitors 
of optically selected QSOs (Sanders et al. 1988). With regard 
to the properties addressed in our comparison study, these PG 
QSOs are indistinguishable from the others and respresent a 
significant contribution to host galaxy studies of radio-quiet PG 
QSOs. 

McLure et al. (1999) conducted an optical imaging study 
of six radio-loud QSOs, nine radio-quiet QSOs and four ra- 
dio galaxies. These objects were selected from a slightly larger 
sample (0. 1 < z < 0.35) previously imaged in the near-IR (Dun- 
lop et al. 1993). Observations were made with HST/WFPC2, 
and targets were centered on the Wide-Field 2 (WF2) chip. The 
F675W filter, which is roughly equivalent to the standard R- 
band, was used. 

Bahcall et al. (1997) pointed HST at the 14 most luminous 
QSOs (z < 0.2) in the Veron-Cetty & Veron (1991) catalogue 
as well as six additional objects (0.2 < z < 0.3) satisfying the 
same luminosity and galactic latitude criteria (M v < -22.9 and 
\b\ > 35°). Observations were centered on the WF3 camera 
of WFPC2 and used filter F606W. In addition to addressing 
host morphology and luminosity, the authors carried out a sim- 
ple analysis of the host environment by identifying all apparent 
companions 7 within 25 projected kpc. Magnitudes and separa- 
tions are presented for 13 such objects. 

As a whole, these five studies contain 38 radio-quiet PG 
QSOs (Table 6). In those cases where an object is included in 
multiple samples, the measured luminosities and morphologies 
are generally in good agreement. There are three ocurrences 
where a host galaxy was called an elliptical by one study and 
given no morphological classification by another. In two cases, 
investigators disagree whether a host is a spiral or an elliptical. 
Thus, between 17 and 19 of the radio-quiet PG QSOs are re- 
ported to lie in elliptical galaxies, 13 to 15 are found in spirals 
and six hosts could not be classified. 

The redshift distribution of this PG sample is comparable to 
that of the 2MASS QSOs. However, observations were made at 
a variety of wavelengths, and differing cosmologies were con- 
sidered. Before comparisons could be made, luminosities first 
had to be translated to a common frame of reference. The /- 
band was chosen because it corresponds closely to the F814W 
filter, minimizing the degree of alteration to the 2MASS QSO 
data. This is especially important for the nuclear components, 
as their SED is not well known. 

Absolute /-band magnitudes were calculated for the various 
samples by first accounting for differing Hq and qo values using 
a form of Equation 1 . If primed values denote the originally 
preferred cosmology, and absolute magnitudes are given in the 

6 As described in Surace, Sanders, & Evans (2001), L,y = L(8-1000/im) and 
Lbbb is the "big blue bump" luminosity at 3200-8400 A. 

7 "Companions" were defined as extended objects with My < -16.5, assuming 
the same redshift as the host galaxy. 
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Table 6 
pg qsos from literature 



Object z Host M Ihost M If/uc Inter- 
(PG) action?* 



N 

N 
N 
N 
Y 

N 
N 
N 
N 
N 
Y 
N 
Y 
Y 
N 
Y 
Y 
N 
N 
Y 
N 
N 



Y 

N 
N 
N 
N 
N 
N 

N 
N 
Y 
N 
Y 
N 
N 

Y 
Y 

N 
N 



* Determined by authors of this paper based on published imaging. The same cri- 
teria applied to the 2MASS sample were used (double-nuclei, merging galaxies 
or disturbed morphologies such as tidal debris and obvious asymmetry). 

a McLeod & McLeod (2001) 

*Percival et al. (2001) 

c Surace, Sanders, & Evans (2001) 

<*McLure et al. (1999) 

f Bahcall et al. (1997) 



X-band, then, 

M x = M x + 5 log[h/h] + 1 .086( 9 o- q' )z. (6) 
The conversion from X-band to /-band is then given by, 

Mi =Mx-2.51og [{ft /fx) {fxfl/fifl)] , (7) 

where fo refers to the flux of an A0V star in the band of interest 
(Bessel 1979; Campins, Rieke & Lebofsky 1985). Values for 
// and fx were measured in the same way as the K-corrections 
described in §3.1. The SDSS composite QSO spectrum only 
extends to the rest /-band and could not be used to transform 
the ZT-band nuclear luminosities in Percival et al. (2001). The 
red-end spectral index of the composite SED was used in this 
case, /„ oc i>~ x - 5 . 

The resulting host and nuclear /-band luminosities are pro- 
vided in Table 6. Of the eight objects included in more than one 
sample, five of them agree within half of a magnitude for both 
the host and the nucleus. Disagreement (< 1.5 magnitudes) 
among the remaining three objects can be understood in part by 
disparate host galaxy morphologies resulting in different fitted 
luminosities. Adopting the first value listed for each object in 
Table 6, the absolute /-band magnitudes of the radio-quiet PG 
host galaxies range from 0.3 L* to 5.2 L* with a mean of 1.5 
L*. The corresponding nuclear values are -21.7 < M/ < -25.4 
and {Mi) =-23.7. 

5. DISCUSSION 

The galaxies in which the 29 2MASS QSOs are found do 
not differ dramatically from those hosting classical PG QSOs. 
Imaging available to date suggests that the latter are both el- 
lipticals and spirals in roughly equal proportion. The 2MASS 
QSO hosts also include both morphological types, though there 
is some evidence that they are spirals more often than not. As 
many as 26% of the radio-quiet PG QSO hosts were not as- 
signed a morphology by at least one investigating team. The 
larger fraction of the 2MASS QSO hosts with ambiguous mor- 
phologies can be explained in part by the greater efforts in- 
vested by some of the comparison studies in modeling the host 
and assigning a morphology. 

The 2MASS QSOs are not found in comparatively underlu- 
minous or unusually bright galaxies. The top panel of Figure 8 
plots nuclear versus host galaxy luminosity for the 2MASS and 
PG QSOs. The horizontal dashed line indicates the absolute I- 
band magnitude of an L* galaxy. The 2MASS QSO host mag- 
nitudes are consistent with the locus formed by the comparison 
sample, and both groups have a mean value of ~1.5 L*. 

On the otherhand, the 2MASS nuclear luminosities are gen- 
erally fainter, ranging over seven magnitudes. Several of the 
brightest objects are as bright as the PG QSOs. The remaining 
majority do not satisfy the traditional criterion for QSO clas- 
sification as indicated by the vertical dashed line in Figure 8. 
Recall from Figure 1 that 2MASS QSOs are all as bright as PG 
QSOs in the K-band. If the 2MASS and PG QSOs also share 
similar intrinsic nuclear /-band luminosities, then between zero 
and seven (mean = 2.7) magnitudes of extinction in the /-band 
are required to shift their luminosities to the average of the PG 
sample, {Mi) =-23.7. 

The bottom panel of Figure 8 plots the 2MASS QSO nuclear 
luminosities versus their Inuc~Ks colors. A Spearman's rank 
correlation test (p = 0.86) confirms the correlation suggested by 
the linear fit (dotted line). Some of the scatter is due to the 
fact that the K s values plotted are total magnitudes including 
both nuclear light and light from the host galaxy. The fact that 
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FIG. 8. — (top) Host galaxy vs nuclear luminosities. Mi for an L* galaxy (-22.3) and the historical Mb < -22.1 + 51og/i QSO criterion are indicated by dashed 
lines, (bottom) Nuclear luminosity versus I-Ks color (Ks includes nuclear and host light). The dotted line is a linear fit to the 2MASS QSO data points only. 



dimmer 2MASS QSO nuclei are also redder is consistent with 
dust obscuration causing the large range in nuclear luminosi- 
ties. There appears to be no relationship between morphology 
and Inuc-Ks or M Inuc , suggesting that any obscuring material 
is local to the nucleus and not dependent on the host galaxy. 
Where suitable color information was available, PG QSO data 
points have been included on the plot. 

Various studies in recent years have suggested that more 
luminous nuclei are preferentially found in elliptical hosts. 
Among the samples considered in this paper, this claim is sup- 
ported by the results of McLure et al. (1999) and contradicted 
by those of Percival et al. (2001). We find no correlation be- 
tween host galaxy morphology and the luminosity of the nu- 
cleus for the 38 radio-quiet PG QSOs, nor for the 2MASS 
QSOs (Xs-band). 

Galaxy interaction, a possible fueling mechanism for AGN, 
is an important parameter to consider when comparing QSO 
samples. However, visual inspection of a host galaxy's dynam- 
ical state is subject to a number of problems. Increased star 
formation triggered by minor and major mergers is more highly 
contrasted in bluer filters, complicating comparisons of stud- 
ies made at different wavelengths. Additionally, the sensitivity 
and resolution of the instrument used for observations affects 
whether or not faint tidal tails are detected or multiple nuclei 
are distinguished. Of course, the subjectivity of assessing de- 
grees of interaction is also sensitive to both the investigator and 
the investigator's working definition of "interaction". 

Tables 2 and 6 list both the 2MASS QSOs and radio-quiet 
PG QSOs as exhibiting or not exhibiting evidence of interac- 
tion. This sometimes uncertain distinction was made by a sin- 
gle individual to remove some degree of subjectivity. Crite- 



ria for interaction include multiple nuclei, obviously merging 
galaxies, remnant debris such as tidal tails, and clear asymme- 
try. Ten of the 2MASS QSOs (34%) were classified as past or 
present interacting systems, as were ten of the radio-quiet PG 
QSOs (29%). Considering the inhomogeneity of the PG imag- 
ing, rough agreement between the two samples seems reason- 
able. 

A more generic indicator of the host galaxy environment is 
the number of companion galaxies within a fixed physical dis- 
tance. In the absence of known redshifts, all apparent com- 
panions (extended objects with small projected separations) can 
be considered for comparative purposes. The methodology de- 
scribed in §3.5 was chosen in order to duplicate the Bahcall 
et al. (1997) study of apparent companions around a sample 
of radio-quiet PG QSOs. Only those objects with projected 
separations of 25 h" 1 kpc or less are presented by those au- 
thors. However, this coincides with the maximum distance at 
which apparent companion information is complete for both 
the 2MASS QSOs and the QDOT ULIRGs drawn from the HST 
archive. Thus, a fair comparison is possible between these three 
samples. 

Figure 9a shows the apparent companion frequencies for 
each sample. The number of 2MASS QSOs, radio-quiet PG 
QSOs and ULIRGs with no apparent companions is 16/29 
(55%), 5/14 (36%), and 1/30 (3%), respectively. Not shown 
on this plot is one ULIRG with five such objects. Table 7 lists 
results of KS tests applied to the frequency distributions as well 
as the projected separations. These values indicate the prob- 
ability (in percentage form) that two samples are drawn from 
the same parent distribution. The number of apparent com- 
panions around the 2MASS QSOs is not statistically different 
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# of companions h" 1 kpc companion separation, r (h 1 kpc) 



FIG. 9. — Comparison of 2MASS QSO, radio-quiet PG QSO and ULIRG apparent companions with M/ < -18.3 and < r < 25 h~ l kpc. (a) Apparent companion 
frequency (fraction of sample with zero, one, two or three projected companions); one ULIRG with five is not shown here, (b) Positions of apparent companions 
relative to QSO or peak brightness. Filled circles correspond to 2MASS QSOs and unfilled circles represent ULIRGs. The large circle indicates the completeness 
limit, r = 25 h' 1 kpc. (c) Cumulative distributions of projected separations for all three samples as well as the expectation for galaxies randomly distributed in the 
plane of the sky. 



than the number around the radio-quiet PG QSOs. However, 
the 2MASS QSO and ULIRG distributions differ at the 99.97% 
confidence level. 

Figure 9b shows the positions of apparent companions rela- 
tive to the central object in both the 2MASS QSO and ULIRG 
samples. The former are indicated by filled circles. Given the 
similar samples sizes (29 and 30), the far greater number of 
apparent companions in the ULIRG images is clear. 

The projected separations are binned and plotted as cumula- 
tive histograms for all three samples in Figure 9c. Additionally, 
the profile expected for objects that are randomly distributed 
in the plane of the sky is shown. Both the 2MASS QSOs and 
the ULIRGs show an excess over random at close separations. 

Table 7 

Kolmogorov-Smirnov statistics* 

Apparent Companion Frequencies 

B97 a ULIRG b 
2MASS C 82.5 0.03 
B97 ••• 0.32 



Apparent Companion Separations 

B97 ULIRG RND d 

2MASS 2.59 50.1 0.01 

B97 ••• 1.60 79.0 

ULIRG 10 7 



* The KS test is a two-sided non-parametric statistic which returns the normally- 
distributed probability that two samples were drawn from the same parent pop- 
ulation. Values in this table are presented as percentages. 

"Bahcall etal. (1997) 

fc QDOT ULIRGs (this paper) 

r 2MASS QSOs (this paper) 

^random distribution 



The cumulative distribution rises sooner for the latter, reflecting 
smaller projected separations of the closest extended objects. 
The probability of being randomly distributed in the plane of 
the sky is 10" 9 and 10" 3 for the ULIRG and the 2MASS QSO 
apparent companions, respectively. The apparent companions 
identified by Bahcall et al. ( 1 997), for their subset of radio-quiet 
PG QSOs, cannot be distinguished from the random expecta- 
tion (with respect to radial distribution only, not total number). 

Interpretation of the comparative apparent companion results 
is not clear-cut. The radial distribution implies a similarity be- 
tween the 2MASS QSO and the ULIRG environments. How- 
ever, if 2MASS QSOs represent evolved ULIRGs, one might 
expect to find more closely separated companions around the 
former. The opposite is in fact seen. The random radial dis- 
tribution of apparent companions around the PG QSOs may or 
may not indicate that previously close companions have already 
been swallowed by the host galaxies. What is clear is that the 
ULIRGs studied have many more apparent companions at all 
separations (within 25 A -1 kpc) than either the 2MASS QSOs 
or the radio-quiet PG QSOs, and unlike the other two samples, 
the ULIRGs nearly universally have at least one apparent com- 
panion. 

The immediate environments of the ULIRGs and 2MASS 
QSOs are significantly different. Farrah et al. (2001) were able 
to classify morphologies for only two of 23 ULIRGs (both ellip- 
ticals). All but one of the objects in their sample are obviously 
interacting systems. Similarly, Surace et al. (1998) found that 
eight of the nine objects in their sample exhibit evidence that 
they are advanced merger systems. These results are consis- 
tent with those of Veilleux, Kim, & Sanders (2002) who found 
that ^73% of their ULIRGs could be adequately fit with a de 
Vaucouleurs profile and all but one of the 118 ULIRGs show 
signs of strong interation. The remains of rapid star formation 
and universal tidal debris that one might expect to see if the 
red QSOs found by 2MASS represent an intermediate stage of 
evolution between ULIRGs and PG QSOs do not appear to be 
present. Clearly, many of the 2MASS QSOs have undergone or 
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are undergoing interaction events. However, many others ap- 
pear to be isolated, undisturbed galaxies. 

There is growing evidence that 2MASS QSOs may simply 
be explained by some form of the proposed unified model for 
AGN. In that scenario, these objects are normal QSOs that are 
heavily obscured along our line of sight. Smith et al. (2002) 
found the sample as a whole to be more highly polarized than 
other QSO samples. We have shown by properly accounting 
for starlight dilution, that the polarization of the nuclear re- 
gion is significantly underestimated by broad-band polarimetry 
in many cases. Furthermore, the magnitude of the corrections 
is correlated to AGN type, with the nucleus to galaxy flux ratio 
generally smaller for narrow-line objects. All of this combined 
with the lack of narrow polarized emission lines in many of the 
highly polarized objects (Smith et al. 2000, 2003; Schmidt et 
al. 2002) is consistent with dust obscuration near the nucleus. 
If obscuring dust were generally located further away, as in a 
galactic disk, one might expect Inuc~Ks to correlate with mor- 
phology. We see no evidence of this. 

The fact that the 2MASS QSOs have similar absolute K s - 
band magnitudes and that their host galaxies do not distinguish 
themselves significantly from those of PG QSOs in the /-band 
implies that the red colors observed are a property of the nu- 
cleus alone. The 2MASS AGN are dominated by type 1 ob- 
jects, and edge-on host galaxies seem to be under-represented. 
Therefore, even the 2MASS AGN sample is probably missing 
extremely obscured objects like those identified by IRAS. In 
these putative highly obscured objects, the AGN may not con- 
tribute significantly in the /-band. In that case, the /-band flux 
will reach a minimum consistent with the host galaxy flux, and 
these objects will not have easily distinguishable J-K$ colors. 
The AGN might become slightly visible in K$, but may still not 
distinguish itself. 

Chandra observations show that many 2MASS AGN are faint 
(probably absorbed) in x-rays (Wilkes et al. 2002), hindering 
the identification of highly obscured AGN in the nearby uni- 
verse. The presence of extinguishing dust very near the nu- 
cleus, as the evidence suggests for the 2MASS QSOs, should 
make these objects bright in the mid-to-far infrared compared 
with their optical/near-IR emission. While this question will 
be examined for a large number of 2MASS AGN with SIRTF, 
identification of many highly obscured QSOs in the local uni- 
verse may have to wait for future, deep all sky surveys at wave- 
lengths longer than two microns. 

6. CONCLUSIONS 

We have observed and analyzed a sample of 29 red QSOs 
selected from the 2MASS near-infrared survey. HST 7WFPC2 
/-band images were used to model the PSF and subtract nuclear 
light to study the host galaxies. We find that 2MASS QSOs live 
in a variety of galaxy types with luminosities ranging between 
0.3 and 4.2 L*. Of the 26 objects with obvious hosts, 9 were 
found in spiral galaxies and 4 in ellipticals. The remaining six- 
teen could not be classified with confidence. It was argued that 
at least half of these may be spirals based on the distribution of 
measured eccentricities. Roughly a third of the sample exhibit 
evidence of interaction, predominantly those with ambiguous 
morphologies. 

The immediate environments of the 2MASS QSOs were in- 
spected by identifying apparent companions with projected sep- 
arations less than 25 h~ l kpc. Less than half the sample has one 
or more apparent companions, with one having two and another 
three. When compared to a sample of radio-quiet PG QSOs, no 



statistically significant difference was found. This is in con- 
trast to the same analysis performed on a sample of ULIRGs 
for which all but one object had at least one apparent compan- 
ion. The projected radial distributions suggest that compan- 
ion galaxies lie closer, on average, to the 2MASS QSOs and 
ULIRGs than the PG QSOs, though data is available for fewer 
objects in the latter sample. 

The 2MASS QSOs were already known to be a relatively 
highly polarized sample. The amount of unpolarized starlight 
from host galaxies was measured and used to estimate polar- 
ization corrections. The resulting increase is substantial and 
depends on AGN type, with type 2s generally having the largest 
corrections. 

Absolute /-band magnitudes of the AGN span the range 
-17.0 to -24.2. A strong correlation between Mj NUC and Inuc~ 
K s color was noted, with redder objects appearing fainter. Nei- 
ther nuclear luminosity nor I NUC - K s was found to correlate 
with morphology. 

No indication was seen that the red 2MASS QSOs represent 
an intermediate evolutionary step between ULIRGs and classi- 
cal QSOs. In light of the similarities between host galaxies and 
environments, the results of this study suggest that the relatively 
faint 2MASS QSO /-band nuclear luminosities are caused by 
obscuration local to the nucleus and that the 2MASS QSOs dif- 
fer from radio-quiet PG QSOs only by their orientation with 
respect to the observer. Thus, the red QSOs being identified 
by 2MASS are likely the missing QSO population predicted to 
exist by x-ray and radio surveys. However, the range of prop- 
erties observed leaves open the possibility that more than one 
phenomenon is responsible for the nature of the large number 
of QSOs discovered by 2MASS. 

We are grateful to Marcia Rieke for providing galaxy SED 
templates that extend into the near-infrared. Thanks to Den- 
nis Zaritsky, Chris Impey, George Rieke, and the referee for in- 
sightful comments and suggestions. This publication makes use 
of data products from the Two Micron All Sky Survey, which 
is a joint project of the University of Massachusetts and the 
Infrared Processing and Analysis Center/California Institute of 
Technology, funded by NASA and the NSF Support for this 
work was provided by NASA through grant number GO-9057 
from the Space Telescope Science Institute, which is operated 
by Association of Universities for Research in Astronomy, In- 
corporated, under NASA contract NAS5-26555. 
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APPENDIX A 
COMMENTS ON INDIVIDUAL OBJECTS 

2M000810+1354-Seemingly isolated spiral galaxy with 
only ~1% of total observed light coming from the nucleus. 

2M005055+2933-Either a modestly inclined barred spiral 
galaxy with visible dust obscuration or an atypical edge-on 
warped disk. Observed light is dominated by host galaxy, and 
the corrected polarization exceeds 13%. 

2M010607+2603-Highly polarized object with very little 
relative light contribution from the morphologically ambiguous 
host galaxy (z ~ 0.4). The surface brightness profile is very 
poorly determined; however, the host must either be a highly 
inclined spiral or an interacting system. The polarization angle 
is closely aligned in either case. 

2M015721 +1 71 2-Asymmetrical host galaxy responsible for 
roughly half of the total light. An extended feature in southeast 
corner is present both before and after PSF-subtraction. A pos- 
sible faint companion or tidal debris is seen within 6 projected 
ft -1 kpc below the object. 

2M022750+7327-Alfhough not visible in Figure 2, an ap- 
parent tidal arm (exceeding 20 h~ l kpc in length) extends from 
the bottom of the host image and arcs to the north. Also, there 
is marginal visual evidence that the host galaxy is in fact two 
merging galaxies of equal size. 

2M023430+2438-Residual image and surface brightness 
profile yield essentially no information about the ^2 L* host 
galaxy at z ~ 0.3. 

2M03245S+774<5-Smoothly distributed emission from host 
galaxy with majority of total observed light originating from 
the nucleus. 

2M034857+7255-Sharesthe distinction (with 2M130700+2338) 
of having reddest near-IR color (J-K s = 3.3) and clearly com- 
prises two merging galaxies. 

2M092049+ 1 903- Apparently isolated spiral galaxy in rela- 
tively highly populated field. Only one apparent companion 
meets the adopted magnitude cutoff, although six additional ex- 
tended objects have projected separations of 25 h~ l kpc or less. 

2M125807+2329-Very faint elongated host galaxy centered 
roughly on the QSO nucleus. 

2M730700+233<S-Elliptical galaxy with a very red near-IR 
color (J-K s = 3.3 for AGN+host). 

2M1 45331 +7353-Clear spiral structure with large observed 
light contribution from host galaxy. 

2M151621+2259-Very faint spiral arms emerging from a 
possible bar. 

2M752757+2257-Disturbedhost showing asymmetry and a 
luminous knot of unknown nature. 

2M754307+7937-Although fairly well described by a r 1 / 4 
profile, residual flux exhibits obvious asymmetry. There is 
a possible alignment of the polarization angle and the fitted 
isophotes, however the degree of polarization is very low. 

2M763700+2227-Host galaxy possesses a disk with multi- 
ple luminous knots within 5 h~ l kpc of nucleus. Not evident in 
Figure 2 is possible tidal debris extending from the host in the 
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northwest direction. The corrected polarization is ^7% and the 
polarization angle is closely aligned with the disk. 

2M163736+2543-An apparently merging system with the 
majority of observed emission from the galaxies. Close inspec- 
tion of the inner region reveals multiple nuclei. 

2M765939+7&34-Innerregion clearly resolved into two dis- 
tinct components. Slightly more than half of the observed lu- 
minosity comes from the nucleus which has a corrected polar- 
ization of ~1 1%. Interestingly, the polarization angle is rotated 
almost exactly 90° from the axis through the two components, 
perhaps suggesting scattering of nuclear light by a close com- 
panion. 

2M770003+27iS-Noresidual flux is visible to reveal the na- 
ture of the ^L* galaxy hosting the highest redshift QSO in the 
sample (z ~ 0.6). The corrected polarization is ^11%. 

2M1 71442 +2602-Grand design spiral with three apparent 
companions. The anomalous light trail is identically present 
in both 400 second integrations. 

2M777559+2(S07-Remarkable asymmetry present in PSF- 
subtracted image. The high polarization and orientation of the 
polarization angle are consistent with the extended feature be- 
ing an ionization cone. 

2M222202+1952-Point-like object used for construction of 
the PSF (this explains the odd appearance of the PSF-subtracted 
image in Figure 2). 

2M222221 +1 95 9-Faint, smoothly distributed residual flux 
with disturbed surface brightness profile. 

2M222554+795S-Inner portion of spiral arms seen in Fig- 
ure 2. Not visible is faint extension of these arms. Only ~1% 
of the total observed light is from the nucleus. 

2M225902+i246-Excellent agreement with r 1//4 profile and 
smooth visual appearance. Recent observations have estab- 
lished that this object's J-K$ color barely falls short of the 
2MASS AGN criterion, though it was included in the full anal- 
ysis. 

2M230304+ 1624-C\&sx\y an interacting system. In addition 
to the extremely close projected separation of the two galaxies, 
a faint luminous tail extends towards the northeast corner of the 
image in Figure 2. 

2M230442+2706-Surface brightness profile described well 
by r 1 / 4 law despite interference beyond 2.5 h~ x kpc by apparent 
companion with ^5 h~ x kpc projected separation. 

2M232745+7624-Faint apparent spiral arms. Inclination 
may have been overestimated if host galaxy is in fact a barred 
spiral. 

2M234449+1221-We,vy little visible residual flux around 
PSF-subtracted nucleus. Well behaved surface brightness pro- 
file is misleading as it probes different objects at different radii. 
The system appears to be interacting. A common envelope 
seems to include the QSO and the extremely close companion 
object. The polarization angle is rotated nearly 90° from the 
axis connecting the two bright spots. 



